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Abstract. 

We have used the ionisation equilibrium equation to derive the electron density in interstellar clouds in the direction to 13 stars. 
A linear relation was found, that allows the determination of the electron density from the Mg I and Mg II column densities 
in diffuse clouds. The comparison of normalised equivalent width of 12 DIBs with the electron density shows that the DIBs 
equivalent width do not change with electron density varying in the range n e = 0.01 4- 2.5 cm~ 3 . Therefore the DIBs carriers (1) 
can be observed only in one ionisation stage, or (2) the DIBs are arising in cloud regions (eg. cores or cloud coronas) for which 
we can not determine the electron density. 
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1. Introduction 

The Diffuse Interstellar Bands (DIBs) are broad absorp- 
tion features, seen in the interstellar medium. There are al- 
most 300 DIBs known in the optical and NIR spectrum 
( |Gala zutdinov et al. (2000)). Despite of over 80 years of in- 
vestigations their nature is still unknown (for a review, see 
|Herb ig (1995)). Many carriers have been proposed as the car- 
riers of DIBs, e.g. solid particles, simple molecules, negative 
atomic ions, carbon chains, fullerens. 

In 1985 Ivan der Zwet & Alla mandola proposed polycyclic 
aromatic hydrocarbons (PAHs) as the source of DIBs. Since 
then it is a popular hypothesis. There are however prob- 
lems with obtaining gas phase laboratory spectra of dehy- 
drogenated and/or ionised PAHs to verify this hypothesis. 



Recently Cox & Spaans (2006) have presented simulations of 
PAH charge state distribution in various environments. In 
clouds with various irradiation and density the fractional abun- 
dances of PAH cations, neutrals and anions changes dramati- 
cally. 

We have used the ionisation equilibrium equation to obtain 
the electron densities in individual clouds. The electron density 
have been compared to the equivalent width of DIBs and the 
CH/CH+ lines. The equivalent width of the CH+ line drops 
with rising n e , but no changes of the DIBs equivalent widths 
are observed. 



Send offprint requests to: P. Gnacihski 



2. Column densities and equivalent widths 

The aim of this paper is to check the dependiences between the 
electron density and the carriers of DIBs. In order to determine 
the electron density we had to measure the column densities of 
two elements in two adjacent ionisation stages. Our target stars 
were stars fulfilling the following criteria: 

- reddened stars of spectral type O or B 

- high resolution Hubble Space Telescope (HST) spectra for 
at least Mgl, Mgll lines are available (both, GHRS and 
STIS spectra were used) 

- hydrogen column densities are available 

- equivalent widths of the chosen DIBs are available 

Column densities of Mg I, Mg II, Si I, Si II, C II, C II* were cal- 
culated from high resolution HST spectra. The spectra from the 
ultraviolet spectral range were downloaded from the HST Data 
Archive. The GHRS spectra taken in the FP-SPLIT mode were 
processed with IRAF tasks poffsets and specalign to achieve 
the final spectrum. The column densities were derived using 
the profile fitting technique. The absorption lines were fitted by 
Voigt profiles. The transitions for which the natural dumping 
constant (T) in not known (Mg II 1240 A doublet, Mg I 1828 
A) were fitted with a Gauss function. The cloud velocities (v), 
Doppler broadening parameters (b) and column densities (N) 
for multiple absorption components were simultaneously fitted 
to the observed spectrum. Both lines of magnesium doublet (at 
1200 A) were also fitted simultaneously - v, b and N were com- 
mon for both lines in the doublet. The wavelengths, oscillator 
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strengths (f) and natural damping constants (T) were adopted 
from |Morton (2003)| 

A convolution with a point spread function (PSF) was 
also performed. The PSF for the GHRS spectrograph con- 
sists of two Gaussian components. The "core" Gaussian has 
a FWHM=1.05 diodes, while the "halo" component has 



rate coefficient (C) (see Table 2) were adopted from 



Shull & van Steenberg (1982) 



FWHM=5.0 diodes ( |Spitzer & Fitzpatrick (1993)> . The rela- 
tive contribution of the "core" and "halo" components into the 
PSF is wavelength dependent and was interpolated from the ta- 
ble given by |Cardellief al. (1990)| The Gaussian PSF for the 
STIS spectrograph depends on wavelength, slit and the mode 
of observations. The tables with FWHM for the combination 
of mode and slit can be found in "STIS Instrument Handbook" 
( |Kim Quijano (2003)| l.The FWHM of the Gaussian PSF was 
wavelength-interpolated from these tables. 

The derived column densities, used to calculate the elec- 
tron density are presented in Table 1. The hydrogen (HI) 



column densities were adopted from Diplas & Savage (1994) 
Molecular hydrogen (H2) column densities come from 



|Rachford et al. (2002)| and Savage et al. (1977) The equivalent 
widths of DIBs and CH/CH+ were kindly supported by Jacek 
Krelowski. 

3. Electron density 

The electron density (n e in cm -3 ) was calculated from the equa- 
tions of ionisation equilibrium for two elements. The first ele- 
ment was Mg, because it is easily observed in two ionisation 
stages. The Mg II column density was determined from the 
1240 A doublet, the Mg I column density was determined from 
the 2026 A, 2852 A or 1827 A line. 

The step rise of dielectron recombination coefficient for 
Mg II with temperature causes the decrease of electron den- 
sity, inferred from Mgl/Mgll, with temperature (Fig. Q]). Such 
behaviour enables calculation of electron density, because the 
curve n e (T e ) from Mgl/Mgll intersects with a curve n e (T e ) from 
another element. The equation of ionisation equilibrium for Mg 
is the following: 



n e N(Mg II) T(Mgn) + n e C{Mg n ) 



(1) 



N(Mg I) a md {Mg 2x ) + a d ie(Mg 2 i) 

Where N(MgII) and N(Mgl) [cm- 2 ] are the column densities 
of ionised and neutral Mg; a ra d{Mg2\) [cm 3 /*] is the radia- 
tive recombination coefficient; a,n e (Mg2i) [crrr'/s] is the di- 
electronic recombination rate; T{Mg\2) [1/s] is the ionisation 
rate of Mgl by UV photons; CiMgn) [cm 3 /s] is the collisional 
ionisation rate. 

Because the coefficients a ra d, adie and C depend from elec- 
tron temperature (T e ) we need an analogous equation for a sec- 
ond element to obtain n e and T e simultaneously. For stars: HD 
24534, HD 203374, HD 206267, HD 209339, HD 210839, the 
second element was Si. The column density of Si I was calcu- 
lated form the 1845 A line, and column density of Si II from 
the 1808 A one. From the intersection of the n e (T) curves from 
Mg and Si (Fig. [T} we have obtained the electron density. 

The F coefficients for Mg and Si were adopted from 
the WJ2 model ( |de Boer et al. (1973)| l. The recombination 
coefficients (a ra d and adie) and the collisional ionisation 



For the stars HD 202904 and HD 160578 the ionisation 
equilibrium was calculated from Mgl/Mgll and CII/CII*. The 
column density of carbon was calculated from the C II 1335 A 
and C II* 1336 A lines, using the profile fitting technique (see 
|Gnaciriski (200 0) for details). The equilibrium between the col- 
lisional excitation and radiative de-excitation of ionised carbon 
is described by the equation (Woo d & Linsky (1997)} : 

N(CIF) n e C{Cn) 



N(CII) a rad (C 2 i) 



(2) 



The radiative de-excitation a r ad(C2i) was adopted from 



Nussbaumer & Storey (1981) The collision rate coeffi- 



cient C(Cn) was adopted from [Wood & Linsky (1997)| and 



Hayes & Nussbauer (1984) 



We have also tryed to use Ca as the second element for ob- 
taining n e . Unfortuatelly, the Cal/Call ionisation equilibrium 
curve (eq. 1) for the stars HD 74455 and HD 149757 does not 
intersect the ionisation equilibrium curve for Mg I/II. The prob- 
lem can be caused by the change of n(CaII)/n(CaI) between 
the edge and the centre of the cloud. |Lepp et al. (198 8) have 
found in their numerical simulations, that n(CaII)/n(CaI) de- 
creases from 4800 at the edge to 160 at the centre of the ( Persei 
cloud. 

4. Results and Discussion 

For a lot of stars Mg is the only element with observations of 
absorption lines of two stages of ionisation. From the column 
densities of neutral and ionised Mg we can only calculate the 
maximal possible electron density nf AX . It is simply the max- 
imum of the n e {T e ) curve. Fortunately this maximum {nf AX ) 
is very well correlated with the electron density n e (Fig. |2). 
The correlation coefficient R=0.89, and the linear correlation 
is following: nf AX = 2.84 • n e . This relation was derived from 
n e = 0.01 -4-2.5 cirT 3 and may not hold for denser or thiner envi- 
ronments. The linear relation nf AX = 2.84 • n e probably reflects 
the fact, that most of the clouds for which we can calculate 
n e have the electron temperature T e ~ 150QK. For stars HD 
24912, HD 74455, HD 91316, HD 141637, HD 147165, HD 
147933, HD 149757 the electron density was calculated using 
the nf AX value and this formula. All derived electron densities 
are presented in Table [3] 

Figure [3] presents the relation between the electron den- 
sity and equivalent widths of the CH and CH+ molecule nor- 
malised to the total hydrogen column density. The equivalent 
widths of CH and CH+ include all doppler components. The 
CH abundance does not change between clouds with various 
electron density. In opposition to CH, the CH+ abundance is 
lower for clouds with large electron densities (more recombi- 
nations). Such behaviour is also illustrated on Fig. @] where the 
theoretical relation between the column densities of Mg I and 
Mg II are presented versus the electron density. 

We have checked that the changes in normalised CH+ 
equivalent width with n e are statistically significant. The points 
on Fig. [3] that are in the direction to the same star (connected 
with a straight line) were replaced by an average n e value from 
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two extreme points. The sample of stars was divided in two 
sets. One with n e < 0.4 cm" 3 and the second one for direc- 
tions with n e > 0.4 cm" 3 . We have calculated the average 
W(CH)/H, of and W(CH+)/H to , and theirs standard deviations 
for stars in both sets. The Student's t-variable was calculated in 
order to check the agreement between the averages for direc- 
tions with n e < 0.4 cm"" 3 and n e > 0.4 cm 3 . The average of 
normalised CH for directions with low and high n e agree with 
significance level 0.7. The average W(CH+)/H, of for directions 
with n e < 0.4 cm -3 and n e > 0.4 cm" 3 differs substantially 
(significance level 0.009). 

Figures [5] and [7] presents equivalent widths of 
DIBs normalised to the total hydrogen column density 
(N(H fof )=N(HI)+2N(H 2 )) ploted versus electron density (n e in 
cm" 3 ). One could expect a drop of DIBs equivalent width with 
n e as seen on Figure [3] for CH+. Unfortunately, none of the 
DIBs bands show a relationship with varying electron density. 
There are two possible explanations for the lack of relation- 
ship between DIBs and n e . First explanation, that the carriers 
of the analysed DIBs may be observed only in one ionisation 
stage. At Figure [4] we can see such behaviour for Mg II. In a 
wide range of observed electron densities (n<,=0.009-2.5 cm" 3 ) 
the column density of Mg II does not change by a considerable 
amount. Such behaviour is also seen for CH on Figure[3] 

The second explanation is that DIBs arise in parts of inter- 
stellar clouds where we observe only one stage of ionisation 
of Mg and other elements. The DIBs can arise in dense cores 
of interstellar clouds, where ionised atoms are hardly observed. 
DIBs may also arise at outer (ionised) parts of the interstellar 
clouds, where neutral elements are absent. For both cases we 
can not calculate the electron density. The hypothesis that DIBs 
carriers are formed in outer regions of interstellar clouds was 
already formulated by Snow & Cohen (1974)[ They observed 
that the 4430A and 5780A DIBs are shallower than expected in 
dense molecular clouds. This result was confirmed by observa- 
tions of the Taurus dark clouds made by Adams onefq/. (19 9 1 ) 
High UV flux on the clouds surface may be responsible for ion- 
ising the DIBs carriers, while the clouds cores are shielded by 
extinction on dust grains. 

5. Conclusions 

The results can be recapitulated as follows: 

1. The electron density in the interstellar clouds was deter- 
mined for 13 lines of sight. 

2. Linear correlation between the electron density and maxi- 
mum possible value of electron density from Mgl/Mgll was 
found nf AX = 2.84 • n e . 

3. The normalised equivalent width of the CH+ line drops 
with rising electron density as expected from the ionisation 
equilibrium. 

4. The normalised equivalent widths of the 5780A, 5797 A, 
6614A, 5850A, 5844A, 6203A, 6270A, 6284A, 6376A, 
6379 A, 6660A, 6 196 A DIBs do not change with electron 
density varying in the range n e = 0.01 4- 2.5 cm" 3 (diffuse 
gas). 
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Table 1. Column densities [cm 2 ] derived from the HST spectra. This column densities were used to calculate the electron 
density (n e ). 
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Table 2. Parameters used in the ionisation equilibrium equation. The ionisation rate T was adopted from the WJ2 model of 
|de Boer et al. (1973) The recombinations rates a rac i and as e as well as the collisional ionisation rate C for Mg, Si and Ca was 
adopted from Shull & van Steenberg (1982)] The a ra d(CII) comes from |Nussbaumer & Storey (1981)| The collision strength 



D.i2(T) was fited by an 8-order polynomial to the data given by Hayes & Nussbauer (1984) 
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Directions with the electron density calculated with the formula 


n e = nf AX /2.84 from the Mg column densities. 
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Table 3. The electron densities derived for target stars. The equivalent widths of the CH + and CH lines were adopted from 
Krelowski et al. (1999) and Krelowski (priv. comm.). These equivalent widths include absorption for all doppler components. 
References: (1)- Rachford7 f al. (2002)| ; (2)- |Savage etal. (1977)| ; (3)- |Diplas & Savage (1994)|; (4)-|Jenkins etal (1986)| ; (5)- 
|Lacouref al. (2005)| ; (6) - N(H 2 ) calculated from the formula : N(H 2 ) = 2.9el9- W(CH); (7) - jBohlin et al. (1978)1 ; • 
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6 ±0.5 


11 ± 1.3 


10 ±2 


96 ± 15 






19 ±2.5 




HD 206267 


89.8 ± 1.6 


222.7 ± 3.6 


44.9 ±2.8 




27.3 ±0.5 


44.2 ± 1.5 


72.9 ±4.1 


199.2 ± 4.6 


25.5 ± 1.5 


36.3 ± 0.9 


117 ± 1.1 


20.3 ± 0.8 


HD 210839 


71.2 ±0.9 


253.1 ±2.7 


60.9 ± 2.4 




30.7 ±1.1 


53.8 ±2.5 


90.2 ±3.3 


482 ± 26 


23.9 ± 2 


55.9 ± 1.3 


147.2 ± 2.5 


25.5 ± 0.9 


HD 24912 


36.1 ±0.6 


200.3 ± 2.4 


29.1 ± 1 


37.1 


20.7 ± 0.4 


22.9 ± 1 


23± 1.1 


197 ±3.5 


12 ± 1.6 


26 ± 1.4 


77.6 ± 2.5 


16 ± 1 


HD 74455 


13 ± 3 


31 ±5 






4.5 ±1 


6.5 ± 1.5 


12 ±3 


105 ± 13 




2.3 ± 1 


17.5 ±2 


1.5 ±0.5 


HD 91316 


17 ± 3 


32 ±5 






4± 1 


15 ±3 


19 ±5 


50 ±8 










HD 141637 


8.1 ±0.8 


78 ±3 






7.3 ± 1 


11 ± 1.5 


11 ±3 


220 ± 14 




3.5 ± 1 


16.5 ±2 




HD 147165 


26.3 ± 4.9 


243.3 ±3.1 


9.9 ± 0.5 




16.5 ±0.5 


18.9 ±0.8 


14± 1.1 


142.6 ±2.1 


9.5 ± 0.5 


20.1 ± 0.3 


60.9 ± 1.2 


8.1 ±0.5 


HD 147933 


50.8 ± 2.4 


208 ± 12.5 


28.1 ±0.8 


20.3 


16.6 ±0.7 


22.8 ± 1.3 


20 ± 1.3 


176.4 ±2.8 


11.6 ±0.8 


25.9 ± 1 


64.6 ± 1.7 


11.8 ±0.9 


HD 149757 


30.5 ± 1.5 


66.4 ± 1.9 


15.7 ± 1.5 


10.8 


1 1 ± 0.5 


14.5 ± 0.8 


13.1 ± 1 


68.2 ± 2 


3.5 ±0.3 


18.7 ±0.5 


40.5 ± 2 


4.2 ± 0.4 



Table 4. Equivalent widths of the diffuse interstellar bands in mA (courtesy of Jacek Krelowski). 




Fig. 1. The n e (T e ) curves from Mgl/Mgll (solid line) and Sil/Sill (dashed line) for the v = -31 km/s cloud in the direction of 
HD 210839. The intersection point of this curves allows us to determine the electron density (n e ). 
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Fig. 3. The equivalent widths of the CH and CH+ spectral lines normalised to the total hydrogen column density are plotted versus 
electron density (n e ). The equivalent width of the CH and CH+ lines includes all intervening doppler components. The solid lines 
connect points from different clouds in the direction of one star. The dashed lines indicates second cloud with undeterminable 
electron density. 
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Fig. 4. Theoretical relation between the magnesium column density and the electron density n e . Solid line represents Mg II 
column density, while the dashed line represents Mg I column density. 
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Fig. 5. The equivalent widths of DIBs normalised to the total hydrogen column density are plotted versus electron density (n e ). 
The solid lines connect points from different clouds in the direction of one star. The dashed lines indicates second cloud with 
undeterminable electron density. 
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Fig. 6. The equivalent widths of DIBs normalised to the total hydrogen column density are plotted versus electron density (n e ). 
The solid lines connect points from different clouds in the direction of one star. The dashed lines indicates second cloud with 
undeterminable electron density. 
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Fig. 7. The equivalent widths of DIBs normalised to the total hydrogen column density are plotted versus electron density (n e ). 
The solid lines connect points from different clouds in the direction of one star. The dashed lines indicates second cloud with 
undeterminable electron density. 



